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Abstract —This paper describes a novel approach for combin¬ 
ing Luby Transform (LT) codes and Network Coding (NC) in 
the context of PowerLine Communications (PLC) smart grid 
networks. Multihop transmissions of LT-encoded data on PLC 
networks are considered and algorithms to combine data at 
relay nodes are proposed. Without the need to decode and then 
re-encode the total received data stream, the relay nodes can 
forward the received data stream while adding at the same time 
their own data. Simulation results are provided confirming the 
good performance of the proposed algorithms. 

Index Tenns —Smart grid, Narrowband PLC, Network coding 
(NC), Fountain codes, LT codes. Multihop communications. 

I. Introduction 

The goal of the smart grid is to associate to the electricity 
transmission and distribution networks communication tech¬ 
nologies, in order to optimize the production, the transport 
and the distribution of the electrical energy. NarrowBand (NB) 
PLC, operating below 500 kHz, is believed to be a natural and 
cost effective infrastructure choice for smart grid communica¬ 
tions HI, as reflected in the adoption of specifications like 
G3-PLC, PRIME, and standards like IEEE 1901.2 and ITU-T 
G.hnem. 

The characteristics of the power lines and their mode 
of operation make the PLC channels non-stationary, highly 
attenuated and very noisy. To ensure reliable transmissions 
on PLC channels, the conventional strategies use a fixed 
rate Eorward Error Correction (EEC) code associated with an 
Automatic Repeat reQuest (ARQ) retransmission mechanism 
in case the EEC code has failed. Since ARQ schemes can lead 
to a high number of retransmissions and acknowledgments, 
erasure codes like fountain codes, which exhibit near-optimal 
overhead, are preferred in replacement of the ARQ mechanism 
El, El, El- Network coding is attractive for the PLC smart 
grid networks with the aim of increasing their performance in 
terms of speed, reliability and efficiency, as demonstrated by 
the extensive literature on this topic El-ii. For a multipoint- 
to-point transmission, performing network coding on several 
LT flows to obtain a larger LT code is not a trivial problem. 
Indeed, LT codes are extremely sensitive to the statistical 
degree properties of the encoded packets. Such properties 
can be altered by a “naive” network coding. The problem of 
distributed LT coding is addressed in this paper for the specific 
NB PLC smart grid architecture. The association of network 
coding and LT coding increases the size of the generated LT 


code and, if optimally combined, the resulting LT code is 
better in terms of overhead. 

The IEEE 34-node power distribution network described in 
El and shown in the Figure [T] reveals that the PLC smart grid 
communication networks presents a tree and radial topology, 
which is considered in this paper. In fact, the network is 



Fig. 1; Node test distribution inspired from the IEEE 34 test 
model. 

modeled as a “line-network” gathering the information from 
the leaves (house meters) to the master (concentrator), as 
presented in the Figure [^(a). 

When polled, each source node sends its data to the con¬ 
centrator either directly or through other intermediate source 
nodes that function as repeaters 0 . A repeating function is 
defined as an integral part of the MAC layer of NB-PLC smart 
grid networks for the purposes of range extension. This type 



(b) 

Fig. 2: (a) PLC smart grid data gathering model, (b) equivalent 
node model presentation as in the IEEE 34 test model 

of network is a “simplified” version of the well known “Y 
networks” configuration, extensively addressed in the literature 
a- Ho). One major difference from the Y-networks, is that 









the relay nodes are forwarding the received LT-encoded data 
stream from a downstream node while adding at the same 
time their own data i.e, all the packets of the relay nodes are 
of degree 1. To the best of our knowledge, this approach has 
barely been addressed in the literature so far. 

The contributions of the paper are the following: 

• The problem of distributed LT codes is analyzed for the 
special architecture of NB PLC smart grid networks, 
where a source node must relay an other source node 
packets while adding its own packets. 

• A relaying algorithm based on the conjoint probability 
on the received degree from Si and the outgoing degree 
from the relay is provided. The performance evaluation 
of the proposed algorithm are presented. 

The rest of this paper is organized as follows: Section II 
presents the system and gives a background on LT codes. Then 
in Section III, we describe the strategies and algorithms of 
network coding at the relays that are retained and we expose 
their design principles. Section IV presents the simulation 
results based on the above schemes confirming the system 
performance. Finally, Section V concludes the paper. 

II. System Model 

Considering the system model of the Figure [T] one can 
notice a repeating pattern in the network. Most of the config¬ 
urations encountered in the NB-PLC networks can be summa¬ 
rized by an elementary “pattern” consisting of one downstream 
source node, a bottleneck relay node and the sink node. Our 
set-up is outlined as follows: 

• The downstream node, called Si, transmits an LT encoded 
message block of length Ki. 

• The upstream node, called S 2 , is collocated with the relay 
( Ri). S 2 is also transmitting a message block of length 
K 2 to the sink. S 2 (Si) and upstream (downstream) node 
will be used interchangeably throughout this paper. 

• When the sink decodes a particular source, it sends an 
acknowledgment to that source. The relay stop relaying 
the discovered packets and the corresponding source stops 
transmitting. 

• There is no buffering for the received packets and there 
is no LT-decoding at the relay nodes. 

In this section a short description of LT codes is given CD. 
Suppose we have to transmit K packets of information. The 
coded packet is produced from the source block M, which 
comprises K source messages M = [mi, m 2 , ■ ■ ■ ,mK], by: 

i) Generating a random variable denoted by dn from the 
predefined degree distribution. 

ii) tn is obtained by the bitwise sum (mod 2) of dn packets 
chosen uniformly at random from the K packets. 

At the decoder side, upon receiving a coded packet, the belief 
propagation LT decoder executes the following algorithm: 

i) If the packet is of degree 1, the packet is considered 
discovered. Then all the previously received and all the 
future packets involving this discovered packet are “xor- 
ed” with the discovered packet. This is to remove its effect 
and to obtain lower degree packets. 


ii) If during the process of step (i) degree 1 packets are 
generated, repeat the step (i). 

iii) repeat step (i) until all the K packets are discovered. 
The performance of the LT decoder is very dependent on 

the degree distribution from which the degrees of the coded 
packets are chosen. This degree distribution is called output 
degree distribution. The optimal output degree distribution 
achieving capacity was found to be the Robust Soliton Distri¬ 
bution (RSD) given by hk'- 


fiKid) 

p{d) 

T{d) 


pjd) + T{d) 


IIK, 

l/{d-{d-l)) 

if d = 1. 

and 

otherwise. 

(2) 

S/K-l/d, 

for d = 1 • • • [AT/S'J — 1. 


S ■ \n{S/6)/K, 

for d = IK/S\. 

(3) 

0, 

otherwise. 



Where K is the number of packets to send, d is the degree 
to be sent, S — c ■ hi{K/S) ■ i/K , the parameters c and <5 
are used to adjust the performance. In the remainder of this 
paper, p.K is a vector of size K + 1 which denotes an RSD 
for K source packets. The indices of px start from zero with 
= 0 . 

For an optimal belief propagation decoding the packets 
must be selected at random with uniform distribution during 
the encoding process. This requirement produces a binomial 
distribution on the number of edges connected to the check 
nodes of the bipartite graph representing the LT code, which 
is called input degree distribution. Although it is easy to satisfy 
this condition for one source, it is quite complicated when the 
source nodes are distributed and network coding is done in 
the intermediate nodes. This important issue is considered in 
this paper. 

For a multihop network consisting of two sources. Si and S 2 
(S 2 being collocated with the relay) with Ki, K 2 information 
packets respectively, the objective is to generate an LT code 
of size K = Ki+ K 2 preserving the input and output degree 
distributions. To respect the input degree distribution at the 
relay node, each time the relay outputs a coded packet of 
degree i, the probability that this packet comprises (i — j) 
packets from S 2 is calculated as follows: 

K2\ 

i<i<K, 0<j<i. 

In order to quantify the impact of a non uniform sam¬ 
pling, some simulations have been carried out. A multipoint- 
to-point transmission is considered with two source nodes, 
K = {1000,200}, Ki = K 2 = Y- The performance in terms 
of the decoding success rate given a specified redundancy 
defined as e = N/K are evaluated, where N is the number of 
packets sent from the relay before the sink is able to decode. 
We assume here that all packets of Si are available at the 
relay. 





At first, all the degree 2, 3 and 4 packets are selected 
exclusively from the source Si or 82 - Therefore a degree 2, 
3 or 4 packet is never coded with the packets coming from 
both sources. The other packets are generated using uniform 
sampling from the overall packets. Then, the degree 5, 6 and 
7 packets are also added to this process and are sampled in 
the same way as the degrees 2, 3 and 4. And then, all the 
packets (of degree less than [^J) are selected either in Si 
or in S 2 . Figure shows the results. The reference curve in 
black, corresponds to the ideal case where we have an LT of 
size K. 

Using this method to generate the coded packets, redundant 
coded packets are more likely to appear. This will induce 
poorer performance of the LT decoding as reflected by the gap 
between the non uniform sampling curves and the reference 
curve. This gap is particularly important, when the number 
of degrees selected in a non uniform manner increases. From 
Figure we note that improvements in terms of overhead can 
be achieved by maintaining a uniform sampling over the set 
of Ki + K 2 = K packets at the relay. 

III. Relaying Strategy 

The objective at the relay is to send packets with an RSD 
of size K while receiving packets with an RSD of size Ki 
and having at its disposal K 2 packets. Furthermore, the relay 
should respect the uniform and random selection of the packets 
resulting in the binomial input degree distribution over all the 
K — Ki + K 2 packets. 

A. Degree distribution management 


This gives a vector of size K +1, obtained as the sum of the 
columns of the matrix P, that is an RSD of size K. 

The marginal probability Psi(<^i = j) is a vector of size 
Ki + 1 that can be obtained as the sum of the rows of the 
matrix P: 

K 

-Psi (j) = X! Flj 0 < j < ATi 

i=0 

Since the received distribution from Si is an RSD, it is not 
always possible to have at the same time a binomial input and 
an RS output degree distributions at the output of the relay. To 
better understand this, let us take a look at the matrix given 
in Q for example. According to this matrix, the relay needs 
to output 22% packets with di = 1 and d = 2. Since we do 
not want to decode at the relay, it is simply not possible to 
have this amount of degree 1 packets from Si because of the 
shape of the RSD at Si. 

We update P to Pq so that the total need for the degree 
j packets coming from Si (PsiO)) is less or equal to the 
available percentage of degree j packets actually coming from 
Si We therefore need to modify the matrix P to 

Po, so that the sum of the columns of Pq (giving Pout) 
must correspond to an RSD and the sum of the rows of 
Po (giving Psi) must always be possible to generate from 
the received Si output degree distribution, in other words 
fsi(j) < MiCiO'), Vj. 

In the following, V(i,fc : n) is an index notation referring 
to the elements Y{i, k),\{i, k + 1), ■ ■ ■ of the matrix 

V. 


Let’s define the matrix P(^_|_i).(ifj^_|_i), with the entry pij 
being the joint probability that the output packet at the relay 
is of degree i and comprises j packets from Si. In an ideal 
case where all the packets are available at the relay the matrix 
P is computed as; 

Pi,j = = i,di = j}. 

= Pr{d = 1 } Pr{di = j \ d = i}. 
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for 0 < j < i. 
for i + 1 < j- 


B. Matrix computation 


Algorithm 1: Compute Pq 


The index begins from zero for notation purposes. For 


follows: 
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0.016 
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4 

0.11 

0.22 

0.11 

0 

0 

(5) 

5 

0.018 

0.058 

0.058 

0.018 

0 


0.005 

0.02 

0.03 

0.02 

0.005 ... 




The relay should therefore send 22% packets of degree 2, 
formed by “xor-ing” one packet from the downstream source 
and one packet from its own packets. 

Since d is the output degree, marginalizing p{d,di) with 
respect to di yields 


Input: 

P{i,j) \> The ideal joint probability of making a degree i 
packet at the relay with j packets coming from Si. 

PKi > The received output degree distribution from Si. 

Output: 

1 * 0 (*ij) A feasible joint probability of making a degree 
i packet at the relay with j packets coming from Si. 

Po = zeros{K + 1, Ki + 1) > Starts Pq to a zero matrix, 
^residual _ ^ Temporary variable for pki- 

^resjduai(o) ^ I ^ Adjust the value of ““"‘(O). 

for i ^ 1 to K do 

Allocate : i) to the elements in Po(i,0 : i) 

proportionally to the probabilities P(i,0 : i) 
maintaining the condition ; 

Ei=r(''^")po{Lj)<FA« 

Update : i) t> Remove the percentage 

allocated for this ith step. 

return Po 


Ki 

F’out(rf = f) = 0<i<K. 

3=0 


To calculate Pq, the elements of the vector pxi dis¬ 
tributed between the different elements of the matrix Pq as 
follows : 
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-All packets are chosen uniformly at random 

” 0.3 

-All packets are chosen uniformly at random 

f -■-■{2,3,4}— degrees packets are chosen in each source 
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Redundancy e = N/K 
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1.2 1.4 1.6 

Redundancy e = N/K 


(a) The number of source packets K = 1000 and K\ — 500 (b) The number of source packets K = 200 and K\ — 100 

Fig. 3: Successful decoding probability of LT codes in terms of overhead when the coded packets are not chosen uniformly at 
random (c = 0.05, d = 0.5) 


As an initialization step, the algorithm starts Pq to a zero 
matrix, hki is put into a temporary variable and 

/r"(0) = 1. 

For the ith row (i starting from 1), the values of Po(i, 0 : 
i) are increased by distributing ^,'^“‘^“'‘*(0 : i) proportion¬ 
ally to the probabilities P(z, 0 : z). This allocation is done 
while respecting the following constraint : 
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Given that it is always possible to generate packets of degree 
di = 0, we arbitrarily set 


residual 

^Ki 


(0) = 1; 


^iraduai reduced, on the indexes 0, 1 , • • • ,i, to 

remove its contribution for this zth stage. 

The main steps of the matrix Pq generation is summarized in 
algorithm 

We present a trace of the execution of the algorithm [T] 
with the following parameters c = 0.05, <5 = 0.5 ana 
Ki = K 2 = 50. The 4 first rows and columns of the matrix Pq 
are computed in each iteration. The 4 first rows and columns 
of the matrix P are computed as follows : 
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^res^iduaij-Q^ being the percentage of packets of degree di = 0, 
coming from Si, employed in the merging process. 

In the first iteration, we allocate : 1) to the 

elements in Po(l,0 : 1). The allocation is done in such a 
way as to best maintain the percentages found in P(1,0 : 1). 
^residual updated on the indexes (0 : 1). 
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In the second iteration, we allocate : 2) to the 

elements in Po(2,0 : 2). As in the previous iteration, the 
allocation is done in such a way as to best maintain the 


Po{i,j) is initialized to a zero matrix and = ^Ki- 



































































percentages found in P(2, 0 : 2). is updated on the 

indexes (0 : 2) afterwards. 
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In the third iteration, we allocate : 3) to the 

elements in Po(3,0 : 3) according to the same process as 
above. is updated accordingly. 
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The update is done the same way for the ith iteration, up 
to i = K. 


C. Output degree generation 

Whenever a degree di — j is received from Si, the relay 
has two options in order to generate the output degree d = i. 

• Either it consults the column with the index j of the ma¬ 
trix Pq and thus forms a pmf (probability mass function) 
by normalizing this column to make it a valid pmf. d is 
drawn from this pmf. It then outputs a degree d packet 
which comprises the received packet and d — j of its own 
source packets sampled at random. 

• Or, it selects to output a coded packet including exclu¬ 
sively its own packets. 

One notes that the relay never receives di — 0 packets. 
Therefore, in order to generate the packets coming exclusively 
from its own set of packets, which happens with probability 
Pexc-si = ^si(O), the relay operates as follows. Upon receiving 
a packet from Si of degree di = j that occurs with the 
probability p-Kiij), at Ps\{j)/p-Ki (j) of the time the received 
packet is mixed with packets of S 2 to form the output packet 
as described before. And for the rest of the time, the column 
of indice zero of the matrix Pq is used to produce packets 
coming exclusively from S 2 . The percentage of the packets 
coming exclusively from S 2 will be the same as Esi(O) while 
the output degree distribution stays RS. 

The detail of the procedure is given in the pseudo code of 
algorithm 

IV. Simulation results 

In this section, simulation results are presented for the 
proposed algorithms and are essentially compared with two 
simulation scenarios; 


Algorithm 2: Packets merging algorithm 

Input: 

Po(d, di) > The joint probability of making a degree d 
coded packet at the relay with di packets coming from 

Si. 

Pcki„ [> The received packet from Si. 

> The received output degree distribution from Si. 

Output: 

Pcko„( > The generated packet to be sent to the sink. 

1 j = degree (Pekin) > The received degree. 

2 T’si(j) = 

3 Uusi„g(j) = Psi(j)/MAi(j) i> The probability of using 
Pekin in the merging process. 

4 if rand < Pusing{j) then 

5 Form a pmf\ 

Pt{D = d) =Po(d,j)/(Ef=oPo('^’i)) for * *e 
choice of the degree d. 

6 Choose d by sampling the pmf. 

7 Sample at random d — j packets from the relay to 
form Pekn. 

8 Send Pekont = Pcki„ 0 Pck^. 

9 else 

10 Choose a degree d according to the pmf: 

Pr{D = d)= Po{d, 0)/( Elo Po{d, 0)). 

11 Sample at random d packets from the relay to form 

_ Pekont. 


• a standard LT i.e. we consider a point-to-point transmis¬ 
sion with only one source that transmits (iTi + K 2 ) LT- 
encoded packets to the sink. 

• a time-multiplexing of two LT codes i.e. the relay alter¬ 
nately sends the LT-encoded packets coming from S 1 and 
its own LT-encoded packets. 

We evaluate the performance in terms of the decoding success 
rate given a specified redundancy defined as e = N/K, where 
N is the number of packets sent from the relay before the sink 
is able to decode. 

Figure shows the performance of the proposed method 
compared with the standard LT and the time-multiplexing 
forwarding scheme for K = 1000 and Ki = 500. It is 
observed that, the proposed merging process outperforms 
the time-multiplexing forwarding scheme, especially when 
the overhead increases. For example, for a 90% successful 
decoding rate, the proposed algorithm requires about 5% less 
overhead than the alternative time-multiplexing LT coding. For 
a smaller size of K {K = 200, Ki = 100), the gap between 
the time-multiplexing approach and the proposed algorithm 
is more important, as described in Figure In this Figure, 
one can see that the curve corresponding to the standard LT 
code is extremely closed to the curve corresponding to the 
proposed relaying algorithm. In addition, the gap between the 
performance of the time multiplexing and the performance 












































Fig. 4: Successful decoding probability of the merging process 
in terms of overhead (with K1 = K2 = 150, c = 0.05,5 = 
0.5) 


of the proposed algorithm is wider; e.g. for 90%successful 
decoding percentage, the proposed algorithm requires about 
10% less overhead than the alternative time-multiplexing LT 
coding. 

V. Conclusion 

This paper presents a strategy of relaying fountain code 
while using inter-session network coding for the special topol¬ 
ogy of NB-PLC networks for smart grid applications. We 
consider the relaying of fountain coded packets on a multihop 
transmission link and we propose algorithms to combine the 
packets at the relay, in a way to preserve the important proper¬ 
ties that allow optimal decoding at the sink. Simulation results 
confirm the good performance of the proposed algorithm for 
a realistic network. 
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